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The Typhoon Ruth and Atmospherics 


By Atsushi KIMPARA 


(Nagoya University) 


Abstract 


The author observed atmospherics concerning the Typhoon Ruth 
which attacked Japan in October, 1951. He found that origins of 
atmospherics are not only scattered over the convergence area of the 
typhoon as in the case of the Typhoon Kezia in September, 1950, but 
also they are distributed in the convergence area in front of the 
trough of westerly waves as well as easterly waves, both of which 
are very likely connected with the typhoon. 

In this case he employed fairly sensitive cathode ray direction 
finders, and he expects to be able to find more detailed characteristics 
of typhoons and other similar meteorological phenomena by using 


more sensitive apparata and waveforms measuring equipments. 


Résumé 


L’auteur a observé les atmosphériques concernant le Typhon 
Ruth qui a attaqué Japon en octobre_1951. Quoiqu’il soit certain que 
les origines des atmosphériques se trouvent dispersés dans la région 
de convergence du typhon de la méme facon au cas du Typhon Kézia 
en septembre 1950, les autres origines se trouvent aussi distribuées 
dans la région de convergence devant la vallée des ondes occidentales. 
dans la haute atmosphére ainsi que celle des ondes orientales, et 
d’ailleurs il est trés semblable que tous les deux ont la relation intime 
avec le typhon concomitant. 
En cette occasion nous avons employé les radio-goniométres a 
rayon cathodique bien sensitives, et nous attendrons que l’on peut 
fournir une méthode nouvelle afin d’étudier les caractéristiques du 
typhon et les autres phénoménes similaires, en employant les appareils 
sensitives pour mesurer la direction ainsi que la forme d’atmosphériques. 
1. Introduction 
- The author derived formerly a general law of generating atmospherics by which 
he explained atmospherics due to the Typhoon Kezia at that time. The characteristics 
of typhoons vary in accordance with seasons of the year and the geographical location. 
In the present article the author intends to describe the principal features of atmos- 
pherics accompanied by the Typhoon Ruth in October, 1951, comparing with those due 
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to the Typhoon Kezia in September, 1950. 
The Typhoon Ruth was found on 6 October, 1951 as a tropical depression on the 


Southern Pacific Ocean east of the Caroline Islands. At 09.00 JMT on 9 October, she 
became a tropical storm and was named “ Ruth” on the sea west of Guam Island; at 
09.00 JMT on 10 October, she developed to typhoon on the sea south of Okino-Tori- 
Shima Island, and reached to the most active period at 21.00 JMT on 12 October, 
indicating 924mb at the centre.“ 

The observation was held from 20.50 JMT on 7 October, 1951 to 20.52 JMT on 
11 October, 1951 at Ohira-Obervatory (35°36/N, 140°28’E) and at Toyokawa-Observatory 
(34°50/N, 137°22’E). The distance between two stations is about 300km and the obser- 
vation was made in accordance with the following programme, taking into account the 
schedule of preparing weather maps: 02.50 to 02.52, 08.50 to 08.52, 14.50 to 14.52, 1750 
to 17.52, and 20.50 to 20.52 JMT. 

At each station we employed cathode ray direction finders with recording 
cameras; 10kc/s was adopted as the observing frequency, taking advantage of its low 
attenuation characteristics in the day-time, smaller night errors, and higher energy 
levels. In order to obtain strictly synchronizing observation, we inserted time marks 
on the film at every second, taken from the standard waves emmitted by the Radio 
Regulatory Agency of Japan. The film speed of every recorder was adjusted to 
10mm/s. 

2. Results of Observation and Interpretation. 

Origins of atmospherics were determined by triangulation based upon observed 
directions of arrival, and we studied the correlation of atmospherics with weather 
phenomena concerning the Typhoon Ruth by comparing with meteorological informa- 
tions in the upper atmosphere as well as on the surface of earth. 

In the case of the Typhoon Kezia we found origins of atmospherics mainly on 
the right hand side of her course, while in the present one, the Typhoon Ruth, we 
found them not only in the places mentioned but also in the converging zone in front 
of the trough of westerly waves in the upper atmosphere such as in 700 or 500mb 
region as well as in the similar region of easterly waves near the centre of the 
typhoon, the latter being, however, not so remarkable as the former. The converging 
region in the trough of easterly waves is generally situated behind it; consequently, 
origins of atmospherics scattered in front of the easterly trough may be one of the 
exceptional cases connected with the typhoon. Moreover, there were lightning flashes, 
showers, cumulo-nimbus. etc., which confirm clearly existence of heavy convergence in 
these districts. 

There is another remarkable distribution of origins of atmospherics off the coast- 
of the Japan Islands on the Pacific Ocean, sometimes prolonged over 180°E, between 
30° and 40°N in a zonal form. This group of origins makes clear itself in autumn 
and winter, while in summer it is quite vague. It seems that this is due to conver- 
gence in the frontal zone in the upper atmosphere as well as on the surface of the 
earth. A jet stream in the upper atmosphere accompanies, in general, polar front 


rl 


The Typhoon Ruth and Atmospherics 38 


under it, and therefore it is very likely that the zonal distribution of origins of atmos- 
pherics resulting from convergence in the frontal zone in the upper atmosphere 
indicates the position of jet stream, but the conclusive description should be postponed 
to future study.* In these observations we found also many groups of origins of 
atmospherics in the southern tropical region, including Philippine and Formosa, cor- 
responding to thunderstorms, showers,cumulo-nimbus, etc., which are, however, inde- 
pendent of the typhoon. Although atmospherics due to cold front in the neighbourhood 
of the Japan Islands are very interesting, they are also independent of the typhoon. 
2-1. Atmospherics due to the trough of easterly waves. 

On 7 or 8 October, when the Ruth was a tropical depression, there was a trough 
of easterly waves in 700mb region which passed through the centre or the west side 
of the Ruth and we observed many origins of atmospherics in front of the trough, i.e. 
on its west side. At 20.50JMT on 11 October, when she developed remarkably as a 
typhoon and indicated 940mb at the centre, there appeared a trough of easterly waves 
in 700mb region in the south-eastern sector of the Ruth, and we observed remarkalbe 
distribution of atmospherics in front’ of the 
trough (Fig. 1). Lightning and cumulo-nimbus 
in the weather map indicate heavy convergence 
there. Generally speaking, there exists a con- 
vergence region in the rear of easterly waves 
and a divergence region in front of it; while 
in summer in 1950 we found very often heavy 
thunderstorms in front of troughs of easterly 


waves. In the Typhoon Ruth we found, more- 


Fig..1 Distribution of origins of atmos- 
over, many origins of atmospherics in front of pherics (x) at 20.50 JMT on Oct. 
LL, 205A. 


troughs in easterly waves accompanying 
thunderstorms, showers, cumulo-nimbus, etc. as an indication of heavy convergence. 
It is, however, not yet clear whether it is due to the typhoon. 

2-2. Atmospherics due to the trough of westerly waves. 

We observed often troughs of westerly waves in 500 or 700mb region which run 
through the centre of the typhoon or near the centre and extend to the north, and 
we found almost always conspicuous distribution of atmospherics in front of the 
troughs. This fact is quite coincident with the general theory that there appears 
convergence in front of the trough of westerly waves and divergence in the rear. 
More over, in our case, the hot and wet air masses accompanied by the typhoon 
entered into the convergence region of the trough and consequently heavy convergence 
phenomena were produced; it is because we found an outstanding distribution of 
atmospherics there. 

For example, at 1450JMT on 9 October (Fig. 2), 08.50 JMT on 10 October, 


* More detailed study will be found in ‘‘ Atmospherics due to Fronts in the upper Atmosphere ’’ 
in this issue. 
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Fig. 2 Distribution of origins of atmos- Fig. 3 Distribution of origins of atmos- 
pherics (xX) at 14.50JMT on pherics (X) at 02.50JMT on 
Oct. 9; 1951. Oct. EY, T95r- 


20.50 JMT on 10 October, and 02.50 JMT on 11 October (Fig. 3), there were troughs of 
westerly waves in 700mb region which run through a line at a distance of 500km in 
front of the centre of the typhoon, and the distribution of atmospherics concerned 
appeared in a region distant about 300km from the centre and scattered along the 
trough. 

At 17.50 JMT on 10 October, the trough 
in 700mb of westerly waves passed through 
near the centre of the typhoon. In this case, 
however, the flow of cold air masses on the 
west side of the trough due to a remarkable 
west wind and the powerful northward flow 
of the wet and warm air masses on the east 
side of the trough due to the typhoon made 
remarkable the distribution of atmospherics 
on the west side and in the neighbourhood 
of the trough, althrough such examples were 

rather rare. (Fig. 4) 

At 0850JMT on 10 October, there 
were troughs of westerly waves in 700mb 
and 500mb regions; the trough on 700mb 


was found far apart in front of the centre, Fig.-4 Distribution Ut exits a seeiaee 
pherics (x) at 17.50 JMT on Oct. 
19% 1951 


and the one in 500mb quite near in front of 


the centre. Both of them extended them- 


yr 


as we 
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Fig. 6 Distribution of origins of atmospher- 
ics (X) at 08.50 JMT on Oct. 18, 1951. 


Fig. 5 Distribution of origins of atmospher- 
ics (x) at 08.50 JMT on Oct. 10, 1951. 


selves far to the north and show a conspicuous distribution of atmospherics in front of 
the troughs in 700mb and 500 mb respectively, suggesting a noticeable convergence over 
there. (Fig. 5). 

There were cases in which we could reasonably presume the existence of the 
troughs in view of the districution of atmospherics as well as isobar, though we did 


not actually find notations of troughs on the weather maps. The following are the 


examples: 
08.50 JMT on 8 Oct., 20.50 JMT on 8 Oct., 08.50 JMT on 9 Oct., 14.50 JMT on 9 


Oct., 17.50 JMT on 9. Oct., 20.50 JMT on 9 Oct., 14.50 JMT on 9 Oct., etc. 


The first two cases of them may be considered to be generated on the conver- 
gence line of wind produced between the centre of the typhoon and the high pressure 
on the north. (Fig. 6). 

In these cases, there were generally cold fronts on the north side of the centre 
and therefore the distribution of atmospherics might be attributed to the convergence 
on cold fronts, but the atmospherics in the convergence area of the troughs of westerly 
waves were too clear to ascribe them to other causes and we could, moreover, found 
thunderstorms, lightning flashes, showers and cumulo-nimbus as marks of heavy 
convergence in these places, whenever we could get meteorological informations there. 
These facts also confirm the law of generating atmospherics published by us 
recently.0@) 

2-3. Atmospherics due to convergence area of the typhoon. 

In the Typhhon Ruth there were also many occasions where we found remarkable 

convergence in 500mb and 700mb regions on the right hand side of the centre, from 


the east sector to the north; some of them may be considered as converging wind 
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lines generated between the centre and the neighbouring high pressure. In these 
convergence areas of the Typhoon Ruth we could find always some outstanding 
distribution of atmospherics almost the same as those in the case of the Typhoon 
Kezia. We can point out easily examples of these kinds in cases of 17.50 JMT on 10 
October, 20.50JMT on 10 October, 08.50 
JMT on 11 October, 14.50 JMT on 11 Octo- 
ber, etc. (Fig. 7). 

3. Conclusion. 

The Typhoon Ruth, which came to 
Japan on the middle of autumn, belongs to 
one of the latest typhoon of the year, and 
consequently her characteristics differ 
somewhat from those of the Typhoon 
Kezia in September. 

As to atmospherics, we observed a 
remarkable distribution of them in the 
convergence area of the typhoon, i.e. mainly 
on the right hand side of the centre as in 
the case of the Typhoon Kezia; we found, 


ics (X) at 08.50 JMT on Oct. 11, 1951. further, another remarkable distribution 

in the convergence area in front of the 

trough of westerly waves in 500mb or 700mb region, and this trough is, moreover, 
often considered to be one of the phenomena to determine the course of the typhoon. 

We found also noticeable distribution of atmospherics in front of the trough of 
easterly waves as well as in the neighbourhood of fronts in the upper atmosphere and 
on the earth in addition to polar fronts concerning jet streams in the upper atmos- 
phere. We shall investigate them in another occosion in connection with atmospherics 
due to cyclones in winter and fronts in the upper atmosphere. 

These observations were made with the instruments of high sensitivity which 
made us enable to investigate the detailed construction of the typhoon through 
atmospherics. We would like to conclude that observation should be made with the 
sensitivity as high as possible so long as it is compatible with the field work, if we 
intend to make radio-meteorology to be useful. 
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Atmospherics due to Fronts in the Upper Atmosphere 


By Atsushi KIMPARA 


(Nagoya University) 


Abstract 


For a long time a distinctive zonal distribution of atmospherics 
has been observed on the Pacific Ocean in lat. 30° to 40°N all the 
year round especially in autumn and winter. Although they are not 
so strong as those from southern districts, they are often found even 
in the interior of Asia at night and have a fairly narrow zonal 
distribution. 

Abundant meteorological informations in the upper atmosphere, 
furnished by airplane observations and radio-soundings in the neigh- 
bourhood of Japan, show that they are generated in the convergence 
region of fronts in the upper atmosphere such as polar fronts. By 
extending this idea to the middle of the Pacific as well as in the 
inland of Asia, where no reliable informations are obtained, we could 
explain with ease the zonal distribution of atmospherics in lat. 30° to 
40°N and long. 80° to 180°EF. 


Résumé 


Pendant long temps une distribution des atmosériques en zone 
distinctive a été observée sur l’Océan Pacifique en lat. 30° a 40°N 
toute l’année, particuliérement en automne et en hiver. Quoiqu’ils ne 
soient pas si forts que les arrivés de la région tropicale, ils se trouvent 
la nuit souvent dans l’intérieur de l’Asie et ont toujours la distribution 
zonal assez étroite. 

Beaucoup des informations météorologiques dans la haute 
atmosphére, obtenues a l’observation par avion ainsi que les radio- 
sondés, indiquent que les atmosphériques précités sont produits dans 
la région de convergence sur les fronts dans la haute atmosphére telles 
que les fronts polaires. En étendant cette idée au centre de l’Océan 
Pacifique et aussi a l’intérieur de l’Asie, ot l’on ne peut trouver riens 
des informations dignes de confiance, nous avons pu commenter sans 
difficulté la distribution zonale des atmosphériques en lat. 30° a 40°N 
et long. 80° a 180°E. 

1. Introduction. 


We have observed for long a remarkable distribution of atmospherics off the 
eastern coast of the Japan Isldnds on the Pacific Ocean in a zonal form in lat. 30° to 
i ee ‘ ; a ivea «4 Geo / ma 
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40°N all the year round. As they are not so strong as those in Idia, the Sunda 
Islands, Australia, etc., we can not find them easily without using sensitive instruments. 
According to a masking effect by strong atmospherics in the neighbourhood of 
Japan it is rather difficult to find these atmospherics in summer; while in spring and 
autumn, especially in winter they appear explicitly in lat. 30° to 40°N and long. 80° 
to 180°E in a zonal form, and in fact they are found mainly in the Pacific Ocean, 
though sometimes specifically at night they are also detected in China, Tibet, Eastern 
Turkistan, etc. ee 
In the neighbourhood of Japan, where we have fairly ample meteorological 
informations in the upper atmosphere, these atmosdherics are found to be scattered 
over the frontal zones in the upper atmosphere and the area between these zones and 
the fronts on the earth. Therefore, if we extend the above idea to regions in the 
middle of the ocean or inland, where any reliable informations are hardly obtained, 
we can understand with ease the existence of these atmospherics by assuming polar 
fronts there, i.e., they are produced in the convergence zone where the ascending 
warm air flow from the south creaps above the cold air masses from the north. 
Really, this kind of atmospherics are 
scattered in a zonal form in lat. 30° to 
40°N, making a long wave which moves 
slowly in the north-south direction, 
revealing certain similarities with 
characteristics of jet streams in the 
upper atmosphere. 
2. Examples and Discussions. 
Fig. 1 shows a distribution of 
Fig. 1 Distribution of origins of atmospherics atmospherics at 20.50 JMT on October, 
Re PE et ets 24, 1951. 11, 1951. A lot of origins of atmos- 


pherics are found to exist along a front in 700mb region, although it is observed as 


a high pressure area in accordance with . 
informations on the earth. A 

Fig. 2 shows a distribution of atmos- 
pherics at 17.50JMT on October 9, 1951. 
Noticeable distribution of origins of atmos- 
pherics are found in the convergence area 
of 700mb region between the front in the 
upper atmosphere and a warm front on the 
earth accompanied by a low pressure. 

Fig. 3 shows a distribution of atmos- 
pherics at 20.50 JMT on October 10, 1951. Re- 
markable distribution of origins of atmos- fig 2 Didiibation of orieits OF aimee. 
pherics are found to be scattered along a pherics (x) at 17.50 JMT on Oct, 
north side of a front in 700mb region, where 9, 1951. 
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Fig. 4 Distsibution of origins of atmospher- 
ics (X) at 02.50 JMT on Oct, 8, 1951, 


Fig. 3 Distribution of origins of atmospher- 
ics (X) at 20.50 JMT on Oct. 10, 1951. 


there is also a low pressure accompanied by cold and warm fronts. 

Fig. 4 shows a distribution of atmospherics at 02.50 JMT on October 8, 1951. 
A number of origins of atmospherics are scattered in the neighbourhood of fronts in 
700mb and 500mb region as well as along a cold front on the earth. 

By observing these examples we may conclude that in the neighbourhood of 
Japan, where we have abundant meteorological informations in the upper atmosphere 
as well as on the earth, the zonal distribution of atmospherics off the coast of Japan 
Islands are located in the convergence region of the frontal zone in the upper atmos- 
phere as well as on the earth. Although in the middle of the Pacific Ocean and the 
inland of Asia there are no reliable meteorological informations, if we extend the 
above reasoning about the distribution of atmospherics, and assume the existence of 
fronts in the upper atmosphere such as polar fronts, we may be able to explain these 
distribution of atmospherics on the Pacific. 

Fig. 5 show a distribution of atmospherics at 02.50 JMT on February 1, 1952. 
Distribution of atmospherics from Kamchatka to the Japan Sea by way of Skhalin 


Fig. 5 Distribution of origins of atmospherics (x) at 02.50 JMT on Feb. 1, 1952, 
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and Hokkaido Island should be attributed mostly to snow showers or snow storms 
there ; while it would be very difficult to explain those on the Pacific scattered in lat. 
35° to 40°N in a zonal form without assuming the existence of convergence in the 
polar front in this district. At 00.00 JMT on Feb. 1 in the neighbourhood of Japan a 
strong westerly air current (110-140 knots) in the upper atmosphere below 500mb passed 
through from lat. 33°-34°N, long. 130°E to lat. 36°-38°N, long. 140°E. The zonal 
distribution of atmospherics coincides fairly well, with the location of westerlies, and 
it seems that the atmospherics exist in the region of strong westerlies exterpolated 


Table 1. 
Passage of ‘‘ westerlies’’ in the neighbourhood of Japan. 


Velocity of wind in the 
Date Time in JMT | Lat.°N in long. 130°E | Lat.°N in long. 140°E region below 500 mb 
(in knot) 
Jan. 28 00.00 35-36 35-36 80- 95 
12.00 35-36 35-36 80- 95 
29 00.00 36-38 ; 36-38 80 
12.00 36-40 38-43 120-145 
30 00.00 36-42 42-44 80-120 
12.00 34-36 36-38 95-125 
Si. 00.00 34-38 36-40 80-125 
12.00 — 3540 36-44 70-100 
Feb. 1 00.00 33-34 | 36-38 80 
12.00 33-35 36-38 — 110-140 
2 00.00 36-38 38-41 100-165 
12.00 34-36 37-40 90-100 
3 00.00 38-40 38-40 85— 95 
12.00 33-35 34-36 100-120 


wid Wea 
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Fig. 6 Movement of zonal distribution of atmospherics from Jan. 28 to Feb. 3, 1952. 


OY es 
an. 28. 08.50 JMT —————————_- X Feb. 1. 02.50 JMT OO - © 

. 29. 14.50 ,, oO Sp Re ce Se SORE RY ey 
SS 80n i450", HH SH Se = = = © 3 oe 06/50. pS ee ee —O 


» 3L 08.50 ,, —--—--—— @ 


12 A, K1IMPARA 


from the above position both on the east and west side, ie. to the continent and to 
the: Pacific. 

-Table I shows the distribution of westerlies near Japan from Jan. 28 to Feb. 3 
with the indication of wind velocities in the region below 500 mb. 

Fig. 6 shows some examples which show very clearly the zonal distribution of 
atmospherics. There we see that from Jan. 28 to 30 the westerlies and the zonal 
distribution of atmospherics move slowly but clearly northward, while after Jan. 30 
both movements do not indicate clear tendency; however, if we inspect in detail, the 
atmospherics returns to the south on Jan. 31 and the westerlies too. On Feb. 1 the 
atmospherics move northward again, and the westerlies too at 12.00 on Jan. 31, and re- 
turn to the south at 00.00 on Feb. 1. On Feb. 2 and 3 the movements of westerlies are quite 
irregular which seems to ascertain the branching and irregular movements of atmos- 
pherics. In fact, the atmospherics moves to the north again on Feb. 1 and on Feb. 2, 
and then it splits into two branches, the one moving to the north the other to the 
south. On Feb. 3 the north branch returns to the south, while the south one stands 
still. All the origins of atmospherics in these days are scattered in lat. 33° to 43°N. 
Taking into consideration of above results, it seems very likely that the zonal distri- 
bution of atmospherics in lat. 30° to 40°N is produced in the convergence region of 
the polar fronts under the jet streams in the westerlies of upper atmosphere. 

4. Conclusion. 

We made extensive observatian of the zonal distribution of atmospherics in lat. 
30° to 40°N and long. 80° to 180°E in its active seasons such as autumn and winter, 
and investigated its correlation with meteorological phenomena by” consulting the 
informations in the upper atmosphere as well as on the earth furnished by the Central 
Meteorological Observatory. 

We found at first the close relation which exists between the distribution of 
atmospherics and the degree of convergence in the upper atmosphere as well as on 
the earth, in studying carefully the informations obtained amply in the neighbourhood 
of Japan; afterwards, we extended cautiously our notion to the case of the Pacific 
Ocean and the inlands of Asia such as Tibet, Eastern Turkistan, Mongolia, etc. where 
the reliable meteorological informations are hardly obtained in general; and referring 
to the results of observation of westerlies in the upper atmosphere near Japan, we 
could finally conclude that this kind of atmospherics is produced in the convergence 
region of the polar fronts, i.e. in the strong ascending warm air flow above the cold 
air masses. The general position of the polar fronts obtained from meteorological 
theory is fairly well coincident with the location of the particular zonal distribution 
of atmospherics above-mentioned. Therefore, further study of position and nature of 
atmospherics will probably reveal us more useful informations concerning behaviour 


and location of the polar fronts and also jet streams related to in the upper atmos- 
phere, | 
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On the Changes in the Atmospheric Electric Field on 
Meteorologically Quiet Days. 


By Minoru KAWANO 


(Electrotechnical Laboratory) 


(Formerly, Geophysical Institute, University of Tokyo) 


Abstract 


The routine observations of the atmospheric electric field, 
atmospheric electric conductivity and space charge were carried out 
from Angust, 1948 to April, 1950 at Hongo, Tokyo. 

In connexion with the electrical characteristics of the compara- 
tively thin layer near the ground, which should play an important 
role for the localization of changes in the atmospheric electric field, 
electric disturbances on meteorologically quiet days are examined. 

The electric disturbances on these days are divided into seven 

"cases, and the correlation coefficients between these electric quantities 

in individual cases are obtained. 

The result shows that, when the time derivative of electric field, 
OE is smaller than 47AE, the electric field are mainly subject to its 
simple relation to conductivity 4, space charge and vertical wind 
velocity w near the ground, namely, AE+pw=const. 

1. Introduction 

It is well known that the changes in the atmospheric electric field on land are 
confined to a rather small area, and this localization seems to be caused by some local 
characteristics in the atmospheric electric conductivity and space charge in a com- 
paratively thin layer near the ground. 

To examine the electrical state of this layer—dust layer—it is convenient to 
discuss the electric disturbances on the meteorologically calm days, because the 
electrical state in the upper atmosphere is considered to change very little in short 
duration on these days, and the interchange of electric charge between the air and 
the earth depends only upon the air-earth current. In this paper, the electrical state 
of this layer is discussed by using the atmospheric electric data obtained at Geophysical 
Institute, University of Tokyo, Hongo, Tokyo. 

2. Instruments 

The routine observations of the atmospher 

carried out continuously from 1947 to 195 


potential gradient, conductivity and spac 


ic electric phenomena at Hongo were 
1, but the simultaneous observations of the 
e charge were carried on only from August 
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1948 to April, 1950. Fig. 1 shows 
the arrngement of the instru- N 
ments used in these observations. 
The observation hut in which 
the instruments were installed 
was shielded from the atmos- 
pheric electric field almost com- 
pletely. 

The instrument used for 
measuring the potential gradient 


consists of a radioactive collector 


and a quadrant electrometer of 


sensitivity of about 10 V/mm/m. ers 

Fig. 1 R: Radioactive Collector 
A: Inhaling Tube of Air 

2.5 meters above the ground, and C: Cage Wire 


The collector was situated at 


its following time was about 
20sec. The reduction factor of this place was 2.2. 

For the measurement of conductivity, a Gerdien-type apparatus having 
capacity of 26cm and a high sensitive Compton type electrometer whose sensitivity 
was about 8mV/mm/m were used. 

The high resistance used was a radioactive resistance of 2x10" and the 
inhaling tube was clinged to the roof of hut, being situated at 2.5 meters on the 
ground. 

The instrument used for the measurement of Space charge was composed of a 
wire cage and a quadrant electrometer. This wire cage which was situated at 1 meter 
high on the ground was 1 meter in its diameter, and a radioactive collector situated 
at the center of the cage was connected to a quadrant electrometer of sensitivity of 
14mV/mm/m. 

In each case of these observations, an optical self-recording device was adopted, 
the time scale of records being 11 mm/hr. 

3. Results of Observation 

Although general feature of the variations of the atmospheric electric field on 
meteorologically calm days does not look much different from the regular diurnal 
variation, remarkable superposing disturbances of short durations and short periods 
were frequently found. 

These phenomena show that there is a definite difference between the meteor- 
ologically calm days and the electrically quiet days. 

Comparing with the simultaneously obtained results of the atmospheric electric 
conductivity and space charge, remarkable disturbances in the atmospheric electric 
field on these days can be divided into the following four cases. 


16 : . M. Kawano 


Table 1 
Electric Field E Space Charge Conductivity 
1 JAE] >| 2] | Ao | large | 4X | large 
2 [AE [> ]£ | | 4p | large | 42. | small 
3. JAE [>{£| | 4e | small | 42] large 
4 [AE [> |£{ | 4p | small | 42 | small 


In this table, | F | indicates the absolute value of the mean of intensity of electric 
field at each time during a whole day, and 4E denotes the deviation of individual 
data from this mean value; do and 4A are the absolute values of deviation from the 
value of mean regular diurnal variation curve. The changes in p and 2 are classified 
either large or small according to their range. Furthermore, the time derivative, 
a3 is also taken into consideration. According as the time derivative of mesa 
field OF is larger or smaller than 47d, classification of remarkable disturbance in 
the electric field in Table 1 can be further divided into the following seven cases. 


Table 2 
hon A Soe | Ap| large | 42] large 
(1-1) | oe | >4m.E (1-2) | ae | <4mE 
ow TRE Se | Ap | large | 42. | small 
(2-1) |r | >40.E (2-2) | oe | <4 
a0 (Arey | Ap | small | 42.| large 
(3-1) | <3 | >40.E (3-2) | | <4n.E 
4. [AE/>]E} | Ao | small [42 [ small 


(4-1) | ae >4m.E 


The examples of these seven cases are given in Fig. 2. 

The relations between EF, o and 4 are given in Fig. 3. In these cases, the values 
of these quantities are the average for 5 minutes. 

The correlation coefficients in each case are given in Table 3. 


Table 3 
y (E, pe) y (E, >) 
1 (1-1) 0.73 —0.71 
(1-2) 0.75 —0.80 
2 (2-1) 0.89 +0.02 
(2-2) 0.86 —~0.01 
3 (3-1) 0.20 —0.66 
(3+2) 0.05 —0.73 
4 (4-1) 0.40 —0.25 


4. Discussion 


If it is assumed that the variations mentioned above is caused by the change of 
electrical state in the dust layer, those variations in the atmospheric electric elements 
mentioned above will be described by the following equations. 
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divE=4z0 (1) 
t=AE+pw (2) 
eS eae 

divi+>,-=0 (3) 


Here w shows the vertical component of wind velocity, and therefore, pw is the 
convection current caused by the air motion near the ground in vertical direction. 
The current flowing in the atmosphere in the vertical direction is the sum of the 
convection current and conduction current. 


From (1), (2) and (3), we get 
1 OF 
div(az4 ow + ini 7) =0 (4) 


As the direction of electric field is generally perpendicular to the ground surface, div 
0 = , : 
may be replaced by Oe" Then, integrating (4), we get 


at 40% 
1E+p0 +-4—Ge=C (5) 


This equation gives the relation between E, p and 4 at every time near the 
ground. When | is sufficiently smaller than 4z74Z, we may neglect the term of 


a i in eq. (5). Then, we obtain the following relations corresponding to each 


case in Table 2. 


1. (2 Mi- E+-AE-2+ 4o-w=0, (6) 
2. (2 4E-i,+4p-w=0, (7) 
3. (2) Ah-E+4E-A=0. (8) 
Modifying these relations we get 
a Re OE 4E 

1. (2) w= — a aes 5) (9) 
. ae A E D 

2. (2) w= oe as (10) 
ual AE Ah 

3. (2) > ie = (11) 


Here, 4E, 4o, 44, 4E ore observable quantities. 

Using the values obtained at Hongo for these quantities; we can estimate the 
values of w by using (9) and (10), where 2# means the vertical component of wind 
velocity U. The horizontal component 7% was directly observed with the Robinson 
anemometer. 


Comparing w thus estimated with observed values of “, we get the correlation 
coefficient in each case as seen in Fig 4 and Table 4. 


Table 4 


1. (2) 1 (w, ”)=0.66 
2. (2) y (w, “#)=0.61 
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Fig. 4 


; These correlation coefficients indicate that w is nearly proportional to u, i.e. 
w=ku. Then, the magnitude of & is estimated as 


1. (2) k=0.02 
2. (2) k=0.03. F 


The coefficients of ratio, — listed above is in a good accord with the values 
usually encountered in meteorology. 

So far as eq. (11) is concerned, the relation between An and = can be as- 
certained directly by electrical observation, and is already given in Table 3, i.e. 
r (E, 4)=—0.73. 

According to these considerations, the relations by (9), (10) and (11) seem to hold 
true. : 

So, we may say that the characteristic and slow changes in the electric field on 
meteorologically calm days can be described sufficiently by eq. (1), (2) and (3), and the 
changes in the electric field are mainly due to the changes of conductivity, space 
charge and wind velocity near the ground. 
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Abstract 


In succession to Part I, the electrical conductivity of the E 
region is calculated in some detail for two atmospheric models. It is 
shown that the factor, by which the tidal oscillation of the E region 
exceeds that at ground level, should be less than 10°, in order to give 
the observed lunar magnetic variation. The calculated lunar vertical 
movement of the FE region is nearly in phase with that observed in 
South England and nearly opposite at Canberra. A possibility is 
shown that the calculated lunar vertical movement of the F2 region 
roughly agrees with that observed. It is suggested that the vertical 
drifts of the F2 region near the magnetic equator may be much 
greater than those in the other latitudes. 
1. Introduction 
In recent years, the electrical conductivity of the ionosphere was discussed by 
many authors. D.F. Martyn [1] suggested in 1948 that the conductivity of the ionos- 
phere might well be enhanced by polarization of the Hall current, as if the earth’s 
magnetic field is cancelled, especially in the D region, and that the D region might 
be the main seat of the currents producing the lunar and solar magnetic variations. 
T.G. Cowling and R. Borger [2] have pointed out that, unless the polarization by Hall 
currents are vertical, this cancellation will be only partial. In two papers [3] and in 
Part I [4], hereafter referred to as I, I have discussed the influence of the Hall 
current concerning the conductivity of the E region, and high conductivity on the 
magnetic equator in the region. I suggested that the E region might be the main 
seat of the S and Z current. Later the problem is discussed by T.G. Cowling [5], K. 
Maeda [6] and D.F. Martyn and W.G. Baker [7] in rather a similar way. The con. 
ductivity and the vertical motion of the ionosphere are examined further in some 
detail in this paper. Two atmospheric models are examined. The height of the 


electro-jet is discussed. The lunar variation is specifically calculated, but some part 


of the conclusions applies also to the dynamo theory of the solar variation. The 


same notation as in I is used, except that time is expressed by ¢ in second and by r+ 
in angular measure i.e. 
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Aca aL, @ = 1.405 x 10-4/sec. 

for the lunar semi-diurnal variation. 
2. General Formulae 

According to M.H. Johnson [8] the diffusion velocity C, in the ionosphere satis- 
fies the equations 

e,n, (E* +e) x H) +e,n,C, x H=(kT/D,)n,C,, (1) 

where r=e, +, —; 
These equations are readily solved for C,, giving the same expressions as Eq. (4) in 
I. When the vertical current vanishes, E, can be expressed by Eq. (15) in J, and, the 


Substitution of this value into Eq. (4) in I gives the vertical component of C,. as 
follows 


y= e-(,— p,- 5)E + (a+ erbe-S )En (2) 
where 

b,=(1/m,»,){cos*$-v,"/(v,2+ w,?)+sin’¢}, 

S,=(1/m,v,)sin ¢ cos $+ ,2/(v,2+,2), 

qr =(1/H) cos ¢-o,?/(v,2+,”), 

P=&)>)n,p,=a'cos*¢+a,sin2¢, 

Q=n,e(q.—4,)=0"'cos 4, 

S=e*>)n,s,=(¢,—0')sin ¢-cos ¢. 


On the other hand, the elimination of E’+c, H, C, and C_ from Eqs. (1) gives 


en,C,=Kp-JI+Kg-J1+Ki-Jxh, (3) 
where | 
Ko=—BiL.L—-1-—NMGtG),  (K.=—Bt-N-L,L.4 MG+L)}.- 
K,=B{l,4+L.4+NLL.—V}, B=L,/{(2,+"L,))\1+N*)}, 
N=L,(L14+4«L2)/(L2+«L), L,=|0;/¥;|=eHD,/kT, Ln=|,/v,|=eHD,/kT, 
«=2A+1. 


Under the assumption that the axes of rotation and magnetism of the earth coincide 
as in I, Ea. (3) gives 
en,W,=K,-cos¢-sin$-J,+ K,-cos$-Jy (4) 


3. The Atmospheric and Ionospheric Models 


The density and temperature distribution in the upper atmosphere has been 
determined by rocket up to 160km height [9]. These values are somewhat different 
from those previously determined by the observations of the meteors and radio wave 
observations, and the reason for the difference is not yet clear. At Pree we cannot 
wholly rely on the rocket measurements, therefore we use the following. two models. 
The mean molecular weight 26 is used throughout. , 

(A) The gas density calculated from rocket observations is per reeted and 
further extrapolated by the formulae of Gledhill [10]. The increase with height of 


the temperature is assumed to be linear, i.e, 
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(5) 


The values of 7 between 120-150km and 150-300km are taken to be 10°K/km and 
6°K/km respectively, the temperature at 120km being taken as 300°K. The density 


distribution is shown by curve A in Fig. 1. 


(B) Using the above mentioned formulae, the model atmosphere by D.R. Bates 
and H.S.W. Massey [11] is approximated. The values of y+ between 90-220km and 
290-300km are taken to be 10°K/km and 6.25°K/km respectively, the temperature at 
120km being taken as 500°K. This density distribution is shown by the curve B in 
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Fig. 1 Atmospheric models (A) and (B). 
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of the same shape, but lowered by 10km 
on the whole, are used for the E and 
F regions. 

To represent the F2 region, we 
take 2,=10*/c.c. in higher 
region than 240km. The above electron 
distributions vs. height are shown in 


constant 


Fig. 2. Using the same diffusion coef- 
ficients as in section 6 in I, we have 
calculated the following values as 


SEC 


Fig. 1. 

In the following, we adopt modified Chap- 
man layers for the scale height with 
linear increase vs. height i.e. A= 
+A(h—h.). For the E region in the 
middle latitudes, the following values are 
used. 

(A) H.o=11km at 4)=120km, and A=0.19, 
(B) Hy=18km at 4o=120km, and A=0.36. 
And for the Fl region in the middle 
latitudes, 

(A) H,.=36km at, 4)>=220km, and A=0.22, 
(B) Hy =54km at 4)=220km, and A=0.36. 
At the equator the electron distributions 


200 


160 


05 a 1.5 Zz 
ELECTRON DENSITY 
Fig. 2 Distribution of electron density with 
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height. Curves A and B refer to 
the models (A) and (B) respectively, 
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functions of height; 
o;'/n;=(e/m,)-v,/(ve+o,?), 
a. /n.=(€/m,)-v./(v¥.2+,"), and a" /n,, 
which are shown in Fig. 3. 
4. The Dynamo Theory 
We treat the case that the region 
between fp, the base of the & region, and 
the upper height /,, makes a uniform 
horizontal tidal oscillation with the velo- 


HEIGHT 


city function 
Vo=ki Pésin(or+a) (6) 


where consideration is specially made for 
o=2 and t=7t’'+A denotes local time in 


angular measure. It was shown by K. 


a es iG” 1 19/9 
Maeda [6] from the equation of continuity ' 4 em.u 
that the order of magnitude of the verti- Fig. 3 Variation with height of o;7/ni, ¢¢7/ne 
cal current density of S, is one hundredth and o!///n, (e.m.u.). Ordinates A and 
: B refer to the models (A) and (B 
of that of the horizontal one. At first we ee oo ®) 
respectively. 


confine our considerations in higher 
latitudes than about-7°. Let E=E’+E’', where E’=c,xH and let E’=(X, Y, Z), then 
we get, using equation ( 19) in I, Z=—cot¢(X+E,,’)—E,‘ =—cot¢-X, which gives under 
the consideration of 7 x E*=0, 

0X/0z=—(Ocot 6/0x)X—cot 6-0X/0x 

0Y/0z=—cot $-0X/ay . 
Since 0B /0x~0E/0y~E/a, it can be shown from the orders of magnitude of (0.X/0z)4h 
and (0Y/0z)4h, that X and Y vary very little for the difference of height ranging from 
100km to 300km. Thus we can use, as a first approximation, the same horizontal 
electrostatic field for the E and F region, in much the same way as in I. 


AY he ee a 
Put a, =\o,dz, B= \ondz, (T=, Ys XY, YX), I= \ Jaz, I= \ Jade 
ho hy ‘ Mo he 


where fh. is taken to be 300km. According to a numerical calculation 8, is almost 


unchanged when /z is higher than 300km. We get the following equations 
T,,=4,{4H,+X)+4(—mH,+ Y), Ih, =B,*X+Bry* Y, 
Ty=4y,(0:H,+X)+a,—mH,+Y), Tey = Byer X+ By: Y. 
where 2,=0¢2/a00, v,=097/asin90A 
The components of J=/,+£ are 
[= 0,°01H,+(4.+ Be) X— Gey H+ (Ary + Bey) Y 


GF) 
p= ayes oH, + (dy+ Byx) X— OywH, + (¢y+ By) Y ase 


i] 
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The elimination of X and Y from these equations and 0X/0A—0(sin@- Y)/06=0, 


‘ a eel Ce gel 0 (sind ) a (See) 
gives ara) — 57h) —o9 00 (Se K, Zy) + OF hie 
0 . 
= (HM, +M,)} +p {Hsin 0 (Mga + Marr) (8) 


where 4=(a,+8,)(Gy+By)+(G2y+ Bry)’ 

K,=4/(ayt+ By), Key=4/(Geyt+Bey), Ky=4/(G.+B8.), Ky.=4/(Oyet+By.)s 

My = (4B xy—ByO2y)/4, Mz = {(ay+ By)a,+(Gry+ Bry) Oxy} /4, 

Ms = (0y(0.+Bx)—Grv(Gye+Bye))/4,  My=(@,8y2—B2Oyx)/4. 
These expressions can be transformed in the same way as in section 3 of I. Then 
the right hand side of the equation (8) is expressed as* 


C ope. cing. dr. @t2)0. ps ot ( 6). P 
ah Tlie hey CL WL P$4,:Ssin(ot+a)+71| c0(0)- PS 


a OP? | 
— 9 cos¢g: 0 ) cos er te (9) 
where ®(0)=3cos?0/{(1+3cos?@)-sin ¢}, 


71= (Bal! —alB™)/{(a' + 8)? +(a + 8), 
re={(a' + Ba’ + (a4 +B )a™}/{(a7 + B24 (al! +p"), 


hy hog 


with a” =\a"dz, pr=\ordz, (r=I, ID. 


ho hy 

For numerical example, two cases are considered, 

Case (i) The model (A) is used with 4,=150km. We get ¥.=0.132-75=0915 

Case (ii) The model (B) is used with #;=150km. We get T1=0:249: 7. =02a7" 
In these cases, the term containing 7, is smaller than one tenth of the term contain- 
ing r2, So that the former can be safely neglected. 
As I is derived from current function R, express R by R= rn Pr sin(mr+on). Then 
the left hand side of equation (8) is reduced to 


i sind. (n+ 1)rnPr sin (mr+ an) + = =r -m-rn Prsin (mr +a™) 
1 cos ¢ m im m. 
ik ak, 1+-3cos20 -2m- *Tn° Prcos(mr+a” ) (10) 
where Ky= (a +82 +(al +R") (a! +), 


Kee {(a? +-A)?-+ (al! + B1)2} (al! + gi), 


For the present cases, the second and third terms of the expression (10) are so small 
compared with the first term that they can be neglected against the first term, leaving 
an error much less than 30%. The equality of the first terms in expressions (9) and 
(10) gives, for o=2, 


d (gms 7°P 3 sin (2r +a) (11) 


" In Eq, (25) of I, the term containing 0Ps/20 is dropped, which cancels mest part of the term 
containing oP (0), ; or Vines 


SGD are 
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where 
wre ‘ 
ri =20-K,-ke, C=—4, 
K,=Ky-re=a' +(a% +B )al™ /(a™ +B"). (12) 


As in I, K, is interpreted as an effective conductivity of the region ranging from ho 
up to 4;. When §”=0, the value is identical with that of (32) in I. And &, is inter- 
preted as effective conductivity of the combined E and F regions with the same tidal 
velocity. We get, for the model (A), A,=3.35x 10°, 

he ho 

\ovdh =8.45 x 107°, \ol'dh =1,45x10-8, (e.m.u.) 

ho ho 
and for the model (B), A,=9.6x 107%, 

ho ho 

\o%dh =5.69 x 10°, \o"'dh=4.72x 10" 

ho ho 
When J, and J, derived from (11), are substituted into Eq. (7), a first approximation 
of X and Y is obtained by 


X =(Ck?/a)(K,/Ki)P 2! {(4/3) sin’¢ —2}cos (27+ a) 
+sin ¢{(2/15)(K./Ko)0P;*/00+cos 0-0P,°/00 -(a!' /K,—a'/K,)}sin (27 +4)] 
=~ (Ck,?/a)-D,-sin(2-+a+ x(X)), (13) 
Y=(Ck,?/a)[2(K./K,)- Pi'-sin (27 +4) 
4+ P,'{(4/3)(K,/K,)sin¢ + (2/sind)(a"/Ki—a"/Ko)}cos(2r + 4)] 


=—(Ch,?/a)- Dysin (2r+0+2(Y))- (14) 
The components of E’=c, = H are 
KM E,' = —(Ch,2/a)3sin20-sin(2t-+a+ jae 
240 - ; 
E,' = —(Ck,"/a)3sin20-cos 6-sin (27+a). 
220 (15) 
5. Results of numerical calculation 
200} The values of a, a!; B, BM vs. hy for 
two models (A) and (B) are plotted in 
180 Fig. 4, under the assumption that 4=0. 
iad The obtained values of K, for various 
160 values of f, and for models (A) and 


(B) are shown in Table 1. The limits 

of errors for K, are estimated to be 
10 o less than 30%. 

Fig. 4 Variation with &, of af, o!7, an d BF, In these ce pe ar be magni- 
err (t= i’ able BF in ie atéh ate) tude of the magnetic field is taken as 
for the models (A) and (B). j H=04 (em.u.). S. Chapman [12] has 

found that to give the observed 


the conductivity integrated through- 


lunar magnetic variations at sunspot maximum, 
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Table 1 


130km , 150 km 180 km 240 km 


oe pak Us? Disk alee 3.0K 10>" 3.3 x 10-8 


1.09 x10 8.1X10? 7.8X 10? 7.6X 10? 


Tinsel =e 2.4x10-8 4.1X10-8 6.8x 10-8 


1.9 «103 1x 103 6.110? 3a K107 


out the effective region must be (2.510-°J-1) (e.m.u.), 2 being the factor by which 
the mean speed of the lunar tide throughout the region exceeds that at ground level. 
The values of 7 are also shown in Table 1. 

The value of AK, should be compared with the value 7x10-°e.m.u. inferred by 
T. Nagata [13] as the integrated conductivity of the ionosphere over Japan from the 
solar flare type variation in geomagnetic field. It is seen that these coincide in the 
order of magnitude. 
6. Conductivity on the magnetic equator 

The conductivity o'!’ vs. height is plotted by full lines in Fig. 5 for the models 
(A) and (B) with maximum electron density at 110km, while dotted line refers to 
maximum electron density at 100km for the model atmosphere (A). In these calcu- 
lations H=0.3 is adopted. 

According to the dynamo KM 


240 
theory, the electromotive force 


in EW direction in these lati- rine 

tudes will be almost electrosta- 200 

tic field. Considering the rela- ae 

tions 0Y/0z=0Z/0y~Z/a it is a 

inferred that Y is almost con- he, 

stant between the E and F @ 49 

regions and the intensity of A ie 3 

is almost proportinal to ¢”” in =~ 

these regions. Referring to the 100 

full lines in Fig. 5, for the 10 

model (A) we get a! 4 QUI 
he Fig. 5 Variation of ¢7Ir vs. height for the models (A) 

=fo"dze=3.4x10-", (al! 4. ptt and (B). 
ho 


/K,=10, and for the model (B) al 4 BUT 710-1, (a! 4 BTT)/K=7.3. Rocket obser- 


vation [14] shows that so-called electro-jet is concentrated between 93km and 105km 
hy he ; 


The ratios jo! ah/\o"" dh, for 4;=120km, calculated by full lines A and B in Fig. 5 
No ho 
are 0.89 and 0.46 respectively. To reconcile with this observation, model (A) is better. 


But further rocket observation is very desirable to confirm this situation, since it has 


succeded only once. D,.R, Bates and H.S.W. Massey -[15], after exhaustive investiga- 
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tion, have shown that A<1 in the E and Fregions. We shall examine the distribution 
of conductivity for the assumed values of 4 The distribution of o7/n, is shown 
with height of the model (A) for «=24+1=1, 5, 10, 20 constant from /) up to 150km 
in Fig 6. Numerical calculation from these KM 


12¢km 150km 
curves gives the ratio of \o"ah/\ ott? dh oe 
: ho ho 

to be 0.91, 0.8, 0.46, 0.41 corresponding ue 
respectively to «=1,5, 10,20. Therefore it 130 
seems appropriate that 4<1. 
7. Distribution of the Electric Field ae 

For numerical example, the follow- D 10 


ing case is considered. Case (iii): The 
model (A) is used with #;=180km. Ac- 100 
cording to Eq. (16) in I, except very near 


the equator, we get a,=a'/sin’¢, a,=0"', 


a,y=a''/sing, and the same form for §,. 
For this case B’%~0, §B/=1.59x 107°, gi= Fig. 6 Variation of o777/n,. with height for 
6.86x10-", of =145x10-". Therefore to  Rivbiciet: as ge meet ee 
get a first approximation, terms containing 

B,’s can be neglected. Further substituting 1 for the denominators sing and sin’¢, 
from Eq. (7) we get for the northern hemisphere 


L,=a'(v,H,+ X)+a'(—mH,+ Y), 
(16) 


I,=—a!"(v,H,+ X)+a'(—mH,+ Y). 
Except near the equator these equations are good approximation, and discussed by 
W.G. Baker and D.F. Martyn [16] for a different ionospheric model. But near the 


equator the discrepancy between Eqs. (16) and (7) is fairly great. From equations 
I,=0R/asin 604, I,= —0R/a00 and (16) and 0X/0A—O(sin0- Y)/00=0, X, Yare eliminated 


giving 
oR afk. 4OR\ _f Ov.) Menten’) | 
oR et oe (sin GG) = aKe\— oy + (OOS (17) 
where K, =a! +(a")*/a" (18) 


Suppose that Eq. (16) is satisfied in the range 0<0<0,. For numerical calculation 
0,=85° is adopted. A general value of the integral of (17) is 


= 2 che KL {P°(0) + Ciqi(9)} sin (27 +4) + C.qi(9) cos(27+a)] (19) 
In the range 7—0:<0<z, in the same way, we get 
R= 2 che, {P3°(0) + Csg2()} sin (27 +a) + Csq2(9) cos (27 +4)] (20) 


where qi(9) and q.(9) are defined in section 4 of I. In the equatorial region 9,<0 
<x—0, we assume that H,=0, then it follows that ¢,,=0, ¢,=9) oy~o'™', Since 
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a,>o"", we take o, as infinite, then it follows that 
X=0, therefore 0S/00=0 (21) 


hy 


For the eastward current we have J,=K,y:Y, where Kp~\o" ‘dh (22) 
ho 


Therefore we get 0R/00=Ky-(0S/sin004). From this equation and (21), S can be 


eliminated giving 
Be ( sind ah) =0 
00\ Kp 00 


Under the assumption that Ky is independent of 0, a general integral of this equation 


(23) 


is obtained as 
R=], Ch*-K,-log( tan 5-) {Cysin (2r-+a)+Cycos(2r-+6)} (24) 
vo 


From Eqs. (16) and (19) we get in the range 0<0<4,. 


X=(Ck,’/a)[{—sin 20+(4/15)(q1/sin 0)(C,+ 7C.)} cos (2r+a) 
+ {(2/15)r -OP3°/00+(4/15)(q1/sin 0)(yC,—C,)} sin (27+a)] 
=—(Ck,?/a)D,, sin (2t+a+ x(X)), (25) 
Y=(Ck,"/a)[{2sin0—(4/15)(q,/sin 0)(C, +7C,)} sin (277+) 
+ {27 sin 20+ (4/15)(q,/sin 0)(7C,—C,)} cos (2r +a)] 
=—(Ck,"/a)D,sin(2t+a+x(Y)), (26) 


ait yar 
where r=a" Jal, 


In the range 7—0,<0<zx we get X and Y in the Same way. In the equatorial region, 
from Eqs. (22) and (24) we get 


Y= —(Chs*/a)(2/15)(K,/Kp)-cosec 0- (Cysin (27+) + C,cos(2r-+a)} 
=—(Ch;*/a)Dysin (20-+0-+ x(Y)). (27) 


From the continuity of the normal component of the current and the tangential 


component of the electric field along 0=0,, z—0,, the constants C,, C, etc. are 
hy 


determined. In this case we get K,=3.75 x 10-3, fo" “dh=3.37107". Since the mag- 
no 


netic field and the electric field are not completely uniform, Ky will be smaller than 
hy 


o! “dh, therefore we take Ky=2K,. For these numerical values we get CQ,=—-—C;= 
No 


0.543, C,= —C,=1.14, C;=—19.76, C,=—10.9. Except very near the equator the terms 
containing C, and C, in Eq. (19) are very small compared with the term containing 
P;*(0), and the current function is expressed almost exactly by the form of the Eq. 
(11). In Figs. 7 (a) and (c), the amplitudes of X, Y and E,’, Ey’ are shown in units 
of —(Ck;*/a) with colatitude, phases x(X), x(Y) are shown in Figs. 7 (b) and (d). For 
cases (i) and (ii), Eqs. (13) and (14) are used, and for the case (iii), Eqs. (25), (26) and (27) 
are used. Except near the equator. cases (i) and (iii) give similar results, while in case 
(ii) polarization field is very small compared with dynamo field FE’, on account of the 
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Fig. 7 Calculated values for cases (i), (ii) and (iii). (a) Amplitudes of EZ,’ and X. 
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Fig. 8 Amplitudes of J, and I, for case (iii) in 
units —(2/15)(Ch,"/a)Ke. 


60° 


8. Vertical Motion of the E Region 
For positive ions, we have 


dissipation in the F region, which has 
The distribu- 
tion of I, and I, with colatitude for case 


considerable conductivity. 


(iii) is shown in Fig. (8), in units of 
—(2/15)(Chs*/a)K,. 


On, /0t=I—a'n,n,—O(n, W,)/02, 


where J=ion production per C.c. per. SEC., 


a! —effective recombination coefficient. The 


terms containing horizontal velocity are neglected, since they are estimated sufficiently 


small. 
have 


On,/Ot=q—aln?—O(n. W,,)/02, 


where g=J/(A+1). 


Under the assumption that 4 is independent of height, using n,=(A+1)n,, we 


(28) 


When vertical current is very small, as in the present case, we have W,~(W,+AW._) 


/(i+A), therefore W,~W,=W for 4<1. 


As’ usual approximate treatment we put 


@n,/0t=0. Let z be the altitude at which ,, determined by Ea. (28) has a particular 


value and let 2’ be the 
the same value. Then, if dz=z—Z 


altitude for which 2,’ determined by Ea. (28) for W,=0, has 
' ig not too large, save close to the maximum 
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electron density, we have 
4z=(0n, W,/0z)/(0q/0z2), (29) 


where the electron distributions for a W,. and for W,=0 are denoted by 2,(z) and 
n,'(2) respectively. Here we neglect the variation of W, with height. When W, is 
not too large, we have 0n,/0z~0n,'/0z and considering n,/~lq/a', from Fa. (29) we 
get 

dz= W,/(2U a! q f(z, ty} (30) 


where q=qof (2, ¢)’, qo being the maximum value of q, and f should not be close to the unity. 
For a Chapman layer--—local scale height 10%cm, a’ =10-'cm’/sec, maximum electron 
density 1.5x10°/c.c-—the influence of W, with semi-diurnal period and amplitude of 
10°cm/sec. is calculated by the method of Kirkpatrick [17]. 

The obtained 4z from the calculated 2,'(z) and n.(z) near noon has been in 
excellent agreement with those obtained by Eq. (30). Therefore Eq. (30) will give 
fairly good estimation for the following examples. When the order of magnitude of 
vertical current is one hundredth of that of horizontal current, it can be shown that 
Eq. (2) gives good approximation to W,. From Eqs. (2) W, can be expressed, in higher 
latitudes than 7° by 


W.,=—cot b-(e/m:){y;/(veit+o,)} (EB, +X)+ (git QAP:i/P.)}(Ey' +Y). (31) 


When 4 is slightly greater than unity W, is markedly affected by the value of 4, but 
W,. is almost unchanged. At lower altitudes than 120km, the second term of the 
right hand Side of this Equation is much smaller than the first term. According to 
S. Chapman [12], to give the lunar magnetic variations for 1905 mean equinox and 
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Fig. 9 (a) Calculated amplitude of W for case (i), 
curves £,, E, E; and F refer to heights 
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mean solstice, @ of Eq. (6) with o=2 should be —17° and —7° respectively and k,?= 
0.3247-a in Eq. (6). Here we take a=—12°. Calculated amplitude of W at various 
heights for case (i) is shown in Fig. 9 with colatitude. /=2% 10? is adopted and FE, 
£, and FE; refer to the heights 120km, 110km and 100km respectively. Substituting 
a’ =10-Scm*/sec, go=200/cm’/sec. and f(z, 4)=0.5 into Eq. (30), we have 


Az=7.08X 10x W, (32) 


by which we may calculate the amplitude P, of dz. The time ¢, of W upward maxi- 
mum after lunar upper transit is shown in Fig. 9 (b) by lunar hours. The limit of 
error of ¢#, is estimated to be about one lunar hour. 

The observed values in South England [18] are P,;=0.93km, f,=11.2 hr. and at 
Canberra [19] P,:=0.19km, #,=5.lhr. The calculated values at 0=40° and at 120km 
level are P.=0.16km, which is about one fourth of that observed in England, and 
4,=2.3 hr. for case (i). For case (ii), the polarization field X is much less than E,' 
and #, is almost determined by the phase of the latter. When ¢, is determined only 
by E,', we get #:=04Ahr. which is independent of latitude. Thus there is a possibility 
that ¢, is nearly in phase with that observed in England. The calculated values at 
6=60° and at 110km for case (i) are P;=0.05km, which is about one fourth of that observ- 
ed at Canberra, and ¢,=0.3hr., which is nearly out of phase with that observed there. 

At present it is doubtful that the recombination coefficient a’ is strictly inde- 
pendent of pressure [20], therefore the above calculation may be too simple to give 
an exact conclusion. Many accurate observations of the height of the EF region in 
middle and low latitudes are very desirable. 

Next, let us codsider using Eq. (4). Since we have 


K,/ Ky —2A-Ly/ + (2A+-1)L,} 


and L,<1 in lower altitudes than 120km, it follows that K,<K, for A<1. Therefore 
the phase of W, is almost determined by that of J, If the flow of lunar current is 
confined in the lower part of the E region, and height gradient of W is neglected, 
from Eq. (4) and (30), it follows that #,=3.4hr. on polar side against the current focus 
and #,=9.4hr. on equatorial side, and these results are nearly opposite phase with 
the observed results. In the present theory it is supposed that the lunar current 
originates in a tidal oscillation in the E region, but the current produced by the 
concomitant polarization electric field flows through all conductive regions i.e. it is 
not confined in the EF region. 

Therefore the above mentioned results and the present theory give different 
conclusions. Since J, is shown in I to be very intense in the FE region on the mag- 
netic equator, therefore according to Eq. (4) W should be very great in the same 
region. This relation will hold for both S and Z current. An accurate measurement 
of the normal E region on the magnetic equator by a modern technique is very 


desirable. 
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9. Vertical Motion of the F2 Region 

In the F2 region W, can be expressed by the second term of the right hand 
side of Eq. (2) as Eq. (40) in I, and we have W,~W,=W. As the life time of elec- 
trons in the F2 region is much longer than that in the E region, in place of Eq. (28), 
we have approximately 

On,/0t = —0(n,W)/0z, (33) 
to discuss semi-diurnal variation. 

Let W=W,cosw(t—t:), W, be independent of height, o=1.405x10-1/sec. ite. 
semi-diurnal (lunar time), then the height of maximum electron density, z,,, may be 
obtained as follows 

Zm=(W,/w) cos {w(t —t2) —2/2} +2p (34) 
To attain the observed amplitude P, of 42, of 2x10!~4%10‘cm, amplitude of W 
should be 28~56cm/sec. 

Calculated amplitude of W and the value of #, for the case (i) are shown by, 
the curves F in Fig. 9 (a) and (b) respectively. Calculated amplitudes are greater 
than those observed by a factor of about 4. This discrepancy and the phase of 4z,, 
are discussed in a later section. 

10. Consideration on the Conductivity of the F region 

H. Alfvén [21] has shown that when an electric field E;, is impressed for an 
incompletely ionized gas in the presence of a transverse magnetic field, the gas asa 
whole is increasingly accelerated in the direction of Ex H in such a way that eo. xH 
cancels EH; and the current decreases exponentially until a stationary state is reached. 
We shall evaluate the role of such a phenomenon in the ionosphere. Let O-xyz be 
taken so that Oz be vertical upward, unit vector in the direction of Oz be k, the 
force pg by the gravity be downward, the density » be a function of height only, and 
both a magnetic field H and an impressed horizontal electric field E,=(X, Y, 0) be 
uniform. Then an equation of motion is obtained by 


a 
0G =I H- (SP + 0g \ie (35) 


Semi-diurnal variation is treated, putting F,=Ey-e, w=1.405x 10-4/sec. © The 
magnitude of E, is given by Eqs. (13) and (14) with Chk,7/a=—43.2 (Z=2 10°). Suppose 
the gas density 2 at z=0 is 2.718x 10°/c.c., temperature T and local scale height H, are 
constant upward from z=0, and do not change when £;, changes. Electron density 
”, (=n,) is 10°/c.c. (const.) in a region 0<z<2H, and zero outside when E,=0, and 
vertical displacement of the charged particles is small compared with H,. At first we 
treat the case when both H and-E; are horizontal ie. H=(H,, 0, 0), F:=(, Y, 0). 
Let ¢o=(u, v, w), then we take u=v=0 and Eq. (35) can be written as 


Ow 0 UT. seni 
eT ad —JvH,—vg—sP (36) . 


Put H,=—H (H>0). and the mean vertical velocity of charged particles Wy referred 
to the rest system O-xyz can be obtained from Eq. (2) as follows 


Piney) ) age ere ee oe 
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W,=W+w=(Y+wH,)/H+w=Y/H (37) 


The equation of continuity becomes 


a (wbe dn “ 
Since in this region we have J,=o’E, where E,=Y+wH, and it may be shown that 
| Y+wH,|<|Y|, therefore substituting oY for J, in Eq. (36), we find the upper limit 
of the variation of p. For the values of the above mentioned Y corresponding to 
lunar current and H=0.3, the order of magnitude of o’ YH, is one thousandth of og, 
therefore the term (@w/0?/) in Ea. (36) can be neglected. For the gas density in the 
F region we have o’=s-p where s=n,-5.09 x 10-'°/H’, then we get o’ YH,=F-p where 
F=Ys-H,. For the present region these relations hold as well. Since F and g are 
supposed to be independent of height, from Eq. (36) under the assumption that, 


P=» (const.) at z=0, we get, p=p, exp {—-M(g+F)z}, 


where M=m/kT. Substituting this value in Eq. (38) we have, as the upper limit of 


WwW, 


w= io are {M(g+F)z+1} +exp {M(g+F)2}1. 


As | F|<g, this becomes approximately 
wH, =H,- {exp (z/H,)—(1+2/H,)} -twFH,/g. 


Put Y=Y,-e, then at z=H, and for H,=100km, we have wH,=5.7 x 10-1Yo-e'e'**/2), 
Thus w and wH, are sufficiently small compared with W and Y respectively and we 
can see that W,~W. Next suppose that the distribution of the electron density and 
the gas density are the same as in the above mentioned case anda uniform magnetic 
field is not horizontal but an impressed electric field is horizontal ie. H=(H,, 0, H,), 
E,=(X, Y, 0). Then the electric current will flow horizontally, affected by the 
secondary polarization field given by Ea. (15) in I 

The equation for w is treated in the same way as Eq. (36) giving the relations 
|wH,| <|Y| and |w| <| WI. In the present gas density, o,, defined in Eq. (16) of J 
can be neglected against o, and o, and the Eas. of motion become 


pa, =0f Y—uH,H. 
a (39) 
oe = —oa,{(X+0H,)H. 
Put H.20,/0=t0 Hioy/0=My X=X cos (t—t,), Y=Y,cosw(t—t), then we have 


Faz stole— —t)+ 242); tan%2=Hy/, 


E,= Y— ZH= Y,: 
V py +o (40) 


FoF sit t:)+2413, tanz:=p,/0. 


E, =X4+v0H,= = Xo° Yuz+w 


In the ionosphere these results will be approximately valid, since both the magnetic 
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field and the electric field are taken as fairly uniform over far greater horizontal 
distance than the thickness of the ionosphere. Therefore these relations are used as 
a first approximation. Let Ox and Oy axes be taken in the south and east directions 
respectively. Then we get, for the vertical drift velocity of charged particles in a 
magnetic field with dip ¢ corresponding to colatitude 8, 

AVS gy ats) (Ho>0) 
= W, cos {w(¢—#.)+ x2} 
where W, = {cos $-¥o/HV1+3co0s°6 } {o/V p72 +07} (41) 


w= Wa 


This consideration may be taken into the E and Fl regions giving the result that w, 
v as well as w are so small that they can be neglected. In the ionosphere, horizontal 
uniformity may be insufficient and the effect of the Coriolis’ force should be taken 
into account. Thus the present treatment will be plausible only as a first approxima- 
tion. 

One of the main current in the F2 region is to dissipate negative (or positive) 
charge accumulated at the current focus in the E region. The air motion in the F2 
region produced by the ampere force will be a whirl shape of the opposite direction 
to the electri: current producing Z variation. In the following we restrict ourselves to 
the case (i). §8,, and By, in section 4 are safely neglected. Divide 8, into two parts 

"= Sordh, B= Sor dh, (r=x, y), then Eq. (7) becomes 


Ls — 0.,U,H, =i (a,+ B,’ )\X— C,yH, toy Y+ B,(X+ vH.), 
L = Oy WH, + Oy .X —OyWH, 3 (Ay + By’) Y+ By’ y= uH.). 


In the middle latitudes the order of magnitude of O2+B,', Gy+By' are Spa"dh+Sy,0'dh 


=7.59 x 10-° (e.m.u.), and those of 8,", By!’ are So'dh= 5.6 x 10-10, Therefore the 
terms containing 8,”, @,'’ are negligibly small compared with other terms. Neglecting 
4 those terms, J,, Z, and X, Y may be determined 
by the above euations with sufficient accuracy. 
08 Thus a first approximation of W may be deter- 
mined by Eqs. (14) and (41). In both cases (i) 
3 ate and (ii) K, will be a little increased, by 
et the decrease of the effective conductivity of 
= 0.4 © 10 
. ~J 
0.2 8 
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the F2 region. For the gas density in the F region we have p,=5.09 x10--n,-(H, 
/H)?/sec. Calculated values of w/Ms1,?+o* with 0 are shown in Fig. 10. 

To obtain a correct value of W we should multiply the value of W obtained in 
section 9 by this factor, and the result agrees fairly well with observed value. We 
have, in place of Eq. (34), 


2m =(W,/w) cos {w(t—te)—2/2+ 72} +20 
=(W,/w) cos w(t—t,')+ 2, (42) 


t,/ is shown in Fig. 11 by lunar hours after lunar upper transit. The Eqs. (33) and 
(34) are valid when lifetime of electrons is sufficiently long and when this is sufficiently 
short, Eq. (28) should be used.. In the latter case 4z,, is nearly in phase with W [17] and 
the phase is earlier by three hours than for Eq. (42). In the F2 region, conditions are 
intermediate of these cases and nearer to Eq. (42). Thus the discrepancy between 
the calculated phase and that observed is about 1-2 hrs., and we can see that they 
are nearly in phase. It seems probable that W on the equator is much greater than 
those in middle latitudes and this relation will be also true for the S current system. 
It is shown by A.A. Weiss [22] that the greater W is, the more the decrease of the 
electron density from the normal Chapman layer in the F2 region in the day time. 
The low density of the F2 region on the magnetic equator [23] may be due to this 
mechanism. 

11. Concluding Remarks 

Conductivity of the E region is calculated for models (A) and (B). To explain 
the direct measurements with rocket, of electro-jet, by the present theory, model (A) 
is better. It seems necessary that the gas density at 100km level is not far from 

10"/c.c. : 

To give the observed lunar magnetic variation by the tidal oscillation of the E 
region below 1km, with A4=0, it is necessary that /=8.1 10" for case (i) and /=10* 
for case (ii). Moreover it is shown to be possible that the effective conductivity of 
the F2 region available for the dissipation of polarization field is smaller than the 
conductivities obtained in section 5 giving the result that / may be smaller than the 
above values. Effective conductivity of the E region for case (i) is about one second 
of that inferred by T. Nagata. Both theoretical estimation by Wilkes [24] and 
observation of wind using radio fading [25] indicate that /=2 x10’. There are no 
great discrepancy between the present estimation and this value. The calculated F2 
region lunar vertical movement seems to be roughly equal to that observed. The 
possibility is shown that the vertical drifts of the £ and F2 regions near the magnetic 
equator are much greater than those in the other latitudes. due to differrent mecha- 


nisms. 
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On the Residual Part of the Geomagnetic S,-Field in the 
Middle and Lower Latitudes during the International 
Polar Year, 1982-33* 


By Hiroshi MAEDA 
(Geophysical Institute, Kyoto University) 


Abstract 


A detailed spherical harmonic analysis of the S,-field shows 
that this field contains not only terms depending on local time but 
fairly conspicuous residual terms. These residual terms seem to 
arise mainly from an effect of the discrepancy between the earth’s 
magnetic and rotational exes upon the dynamo-action in the S,-layer. 

By coordinate transformation to make the residuals least, we 
get a more suitable coordinate for the S,-field than the geomagnetic 


one. 


1. Introduction 
On the statistical investigations of S, based on data during the Polar Year, 


1932-33, N.P. Benkova [1] has given the variations of S, with longitude (for the 
summer season) being expressed analytically in terms of spherical harmonics. M. 
Hasegawa and M. Ota [2] have investigated the whole data in more detail, and 
indicated the S, current-system at every successive two’ (universal) hours. Accordingly, 
they expressed this changing S,-field in a general analytical formula containing the 
coefficients which vary with both times and longitude, and this field can be divided 


into two parts, the one depending on local time and the other on universal time 


(residual). 

That these residual terms arise from the o 
with respect to its rotational axis, on the basis of the danam 
out by A. Schuster [3], and many workers [4, 5, 6] have since discussed this problem. 
Recently T. Nagata and M. Sugiura [7] have calculated the S,-field taking account of 
the discrepancy between above two axes, based on the dynamo-theory, and obtained 
a result approximating the Hasegawa and Ota’s figures (2}. . 

Because of these residual terms, it seems that the geo 
the entire Earth proposed by A.G. McNish [8] is not always suitable for the S,-field. 

In the present paper, the author has described a detailed feature of the residual 
part of the S,-field, and tried, by calculation to obtain a coordinate system for the 
S,-field which would most probably reduce the residual terms as a whole. 

The data used in the present analysis, which was prepared and analyzed by 


pliquity of the earth’s magnetic axis 
o-theory, was first pointed 


magnetic coordinate for 


o Prof, M, Hasegawa on his sixtieth birthday, 
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Table 1—List of magnetic observatories in the geomagnetic middle and lower latitudes. 


Observatory Abb. ¢g 3 ? A v o At 
Sitka Si 57.0°) 224.7°| 60.0°| 275.4°] 21.4°} 59.9°| 287.7° 
Eskdalemuir Es 55.3 |. 85628 58.5 82.9 | —20.4 Byer 81.7 
Lovo Lo 59.4 17.8 58.1} 105.8 | —22.1 58.0 1m 104.1 
Sloutzk Sl 59.7 30.5 56.0 | 117.5} —20.6 Bisel al6.0 
Rude Skov RS 55.8 12.4 55.8 98.5 | —20.6 55.6 97.4 
Agincourt Ag 43.8 | 280.7 55.0 |_ 347.0 3.6 51-2) Wase9 
Abinger Ab Sls2 I) 359.6 54.0 83.3 | —18.4 52.6 83.0 
De Bilt DB 52.1 Lifer 53.8 89.6 | —18.9 52.8 88.9 
Manhay Ma 50.3 5.7 5200 88.8 | —18.2 ren lea! 88.9 
Val Joyeux VJ 48.8 A i Wes 84.5 | —17.5 50.2 84.7 
Swider Sw OZ at Zee 50.6 | 104.6 } —18.3 51.0 | 104.6 
Nantes Na 47.2 | 358.4 50.5 80.1 | —16.8 49.0 80.6 
Cheltenham Ch SO selighee aed, 50:1.) 350.'5 2.4 46.2 | 356.8 
Kasan Kn 55.8 49.1 49.2! 130.6} —15.6 51.2) “1st 
Vissokaya Doubrawal VD 56.7 60.1 46.5] 140.7) —13"3 51.0 | 141.9 
Vienna Vi 48.2 16.2 47.9 822 ST722 47.6 98.3 
O’Gyalla OG 47.9 18.2 47.1 99.8 | —17.0 47.1} 100.3 
San Miguel SM 37.8 | 334.4 45.6 AG i es ew 42.5 D3..3: | 
Ebro Eb 40.8 0.5 43.9 ty —1500 42.4 Sits 
Zouy Zo 52.50.1040 41.0 "re B74 Ao ag 45.04: 1783384 = 
San Fernando SF ceraeisia) (mesa re ie: 41.0 $1.2) eG 38.8 01 
Tucson Tu 5222) 249.2 A074" 312.2 10.1 38.1] 319.9 
Toyohara Ty 47.0 | 142.8 36.9 | 203.5 G.7 40.8 | 212.7 
San Juan SJ 18.4 | 293.9 29.9 3.2. ies 25.8 O23 ee 
Teoloyucan Te 19.8 | 260.8 29.6 | 327.0 6.6 26.6} 333.4 
Helwan He 29.9 rot es) 2?7.2| 106.4 |" 49°97 Zr. 7 | 1Oor SS ee 
Kakioka Ka 36127)" 140.2 26.0 | 206.0 6.2 29.9 | 211.5 
Aso As 5249.)) 2310 22.0 | 198.0 4:2 26.0] 203.3 
Honolulu Ho 27.30) 20ds9 21.1] 266.5 12.3 21.5 | 27229 
Dehra Dun DD 30.3 78.0 20.5 | 149.9 "6 23.0 peasy aoe 
Lukiapang Lu SS") e170 2030) SOe 7 2.1 24.0 | 194.0 
Z6-Sé ZS po Kees Wea eh IB 19.8 | 189.2 2.2 23.8 | 194.2 
Alibag Al 18.6 72.9 9.5] 143.6] —~ 7.9 2.27 438.0T 
Antipolo At 14.6 | 121.2 3.3] 189.8 2.0 7.3 | 194.8 
Huancayo Ha |) s=8200| 28t7 ) den Wehees LS. = 25 ooeey, 
Ellisabethville El oat BINS 27.5 | =12.8 94.1 | —11.8 | —13.1 IBM Fn foe 
Apia Ap | 13.8) 188.2 | —16.0} 260.2} 11.7] —14.7| 264.1 
Batavia Ba j — 6.2] 106.8] —17,.6| 175.6] — 0.9 —13.7] 180.8 
Pilar Pi S17) 296.4 | Sons 4.6} — 1.1] —24.3 9.4 hes 
Tananarivo Tn | —18.9 47.5 | —23.7]. 112.4 | 11.9 227 109) ae 
Mauritius Mu | —20.1 57.6 | —26.6] 122.4] —10.3 —24.9) 129.2) _ 
Cape Town CT | -88.9) 18.5 | =99.7 | 79.9 | 38.7 agg el] soeme hun 
Watheroo Wa | —30.3) 115.9] —41.8] 185.6 1.3 | —37.7 | 190.8 
Toolangi To | —37.5] 145.5 | —46.7] 220.8 9.5 | —43.1 | 223.7 
Christchurch Cr | —43.5] 172.6 | —48.8| 252.6 15.2 | —46.1] 253.6 
Orcadas del Sud OS | —60.8] 315.0 | 50.09 18.00). 269 395 eae 24; 00 


Notes: 4, }, =Geagraphic latitude and longitude, 
Y =Angular difference in direction at ob 


observatory between the geographic and o 

Te ei q 

meridians, Positive when measured from north around by east “jadi hig t 
: Sarees Pe ep aR 


ae 
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Prof. M. Hasegawa and Dr. M. Ota, are the mean state of summer and winter (say 
¥4(S+W)) of S; as recorded at 46 stations (Table 1 and Fig. 1) in the geomagnetic 


thiddle and lower latitudes (between 60°N and S) during the Second International Polar 
Year, 1932-33. 
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Fig. 1—Distribution of magnetic observatories in the middle and lower latitudes. 
Geomagnetic coordinate. 


2. The residual part of the S,-Field 
As the mean state from the geomagnetic X- and Y-components, a magnetic 
potential V of S, at the earth’s surface can be expressed in terms of the series [2] 


V=r Od, AvP” (cos 9) cos (1T + mA) + ,BrP% (cos 9) sin (IT+ mA) 
nmel 
+,C™P™ (cos @) cos (LT —mA)+,D%P® (cos 9)sin({T—mA)), (1) 


where 7)=radius of the Earth, 
@=geomagnetic colatitude, 
A=geomagnetic east longitude, 
T=—69° meridian time (say geomagnetic universal time, and also T+ A=t 
is geomagnetic local time), 
P™(cos@)=Schmidt’s normalized spherical functions. 


In the series (1), the terms of /=m in the terms ,A7Pxcos(/T+ mA) and ,BuPrsin?T 
+mA) vary with local time (say local time field), and the others do not (say residual field) 
in which the terms of »=0 depend on universal time only (say pure universal time field). 

Tables 2A and 2B give the values of the coefficients ,A7, ,Bn, ,Cn and ,D% in 
(1), in units of force (10-°C.G.S.). The figures in Table 2A show a slight difference 
from those in [2], because the former corrected by recalculations in a different way, 


From Table 2A it is interesting to note that the six terms, 
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Table 2A—Values of coefficients ,A™ and ,B™ in the series (1) expressed in units of force 


H. MaErpa 


(10-°CGS=1y), where ( ) indicates main six residual terms. 


Pe Pe | Pe Pe | Ps | Pe | 
jAo}) -—0.13 0.37 0.84,| —0.12 eee 0.37 | cosT 
oB2 0.10 0.15 wer —0.04 0.66 0.16 | sinT 
2AQ 0.02 0.16 0.25 0.09 | —0.17 | -—0.41 | cos2T 
2Be 0.02 —0.03 0.13 —0.06 0.33 —0.24 sin2T 
P P} P: P: P: Py 
143 | —0.14 8.40 | —0.63 | —1.75 0.60 0.21 | cos(T+A) 
1B} 0.05 | ~—2.35 | —0.45 0.40 0.24 0.12 | sin(T+A) 
2A (ny —0.06 oat 0.14 0.23 0.07 | cos(2T-+ A) 
2B: | \—0.37 0.21 0.69 0.23 | —0.34 | —0.19 | sin@7+A) 
3A 0.13 —0.05 0.01 —0.02 —0.07 —0.01 cos(3T+ A) 
3Bh 0.14;]' 0,10 | “soca 0.08 | —0.04 0.04 | sin(83T+ A) 
P2 P2 P2 P2 P2 P2 
1A2 ( a] 0.06 0.28 0.08 | —0.06 | —0.11 | cos(7+2A) 
1B3 0.91/| —0.41 0.02 | —0.06 0.20 | —0.11 | sin(T+2A) 
yAz 0.00 | —4.15 0.28 0.14 | —0.07 | —0.02 | cos@7+2A) 
Bt 0.76 1.69 0:81}? 20.48 0.09 0.22 | sin(2T42A) 
s4z| 0.00] 0.04] 0.26] —0.0¢] 0.04] 0.02 | cosT+2A) 
3B3 0.27) “=0.18 | =0.80°7 20:30 °F oetiee 0.08 | sin(@7T+2A) 
P3 Ps P3 P3 P3 Ps 
1A3 0.19 0.28 | —0.17 0.14 —— —— | cos(T+S3A) 
1B3 0.26 | —0.15 0.05 | —0.11 — — | sin(T+34A) 
2Ag | = —0.08 | —0.11] -0.01| —] 2°) Geereaay 
Bi —0.61 0.21 —0.15 0.12 —— — sin (2T+3A) 
a are ares eear sy rs Gee ree rere ese 
3Ba | —0.54 | —0.86 —0.04 | —0.02 — —— | sin(8T7+3A) 


Table 2B—Values of coefficients ,C” and ,D™ in the 


series (1) expressed in units of force 


(10-5CGS=1y). 
ae ts ak th i a | Ps 
1Ch} 0.27 | 0.20] -o.17| 0.12 | —0,01 | —0.12 | cos(T—A) 
iDi | —0.43 0.02 0.07 | —0.18 | —0.04 0.03 | sin(T—A) 
2Ci| —0.08 | -0.02 | 0.01 | —0,06 | -0.12| 0.05 | cos(2T—A) 
Di} 0.26 | 0.06 | -0.09 | —0.01] 0.07 | —0.06 | sin@7—A) 
sCi| —0.04 | —0.01| -0.05| — | —J| 0.02 | cos@7T—A) 
1 -€:00'| 0.01 | -0,02 degen 


2D} 


sin (3T— A) 
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cos(T—2A) 
sin Am 2 ee 


cos(2T—2A) 
sin (2T—2A) 


cos(3T--2A) 
sin (8T—2A) 


cos(T—3A) 
sin(T—3A) 


cos(2T—3 A) 
sin (2T—3A) 


cos(3T—3A) 
sin (3T—3A) 


Table 3—Calculated values (by dynamo-theory) of a 7,A™ and ,B™ expressed in 
units of force (10-*CGS=1y). 


cos(T+ A) 
sin(T+ A) 


 eos(2T + A) A) I 
sin (T+ A) 


sin (37+ A) 


cos(T+2A) 
eek ude ia 


cos(2T +24) 
sin(2T+2A) 


cos(3T+2A) 
sin(3T+2A) 


cos(T+3A) 
sin (T+3A) 
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2As | ;—0.90 cos(2T+3A) 
2B ( eats sin(2T+3A) 
oe ee eee tee ar peoreae ye 
3By —1.11 sin (37+3A) 


Table 4—Values of amplitude 2a%={(,A™)2 4 (,B™)2}2 expressed in units of force (10-°CGS 
=1y), and contribution of velocity potential to magnetic potential in the dynamo- 
theory, where ( ) indicates main six residual terms. 


ai 0.16 0.40 ARP 0.13 0.91 \| 0.40 i re 
a 3 Car 3,0} (yiu2) gig3 | Sines 
go} 9-03 | 0.16 | 0.28 | 0.11 | 0.87 | 0.26 | _ ' 
2 2 | v8 | vans | van2 | vaya | eT +22) 
Pid Pi} ok.) Pas te Pt ae 
oe OFS 8.74 Or? 1.79 0.65 0.24 c 
"| ga | oben | ysgs | oten | uaa | vags | Sim+ A tia 
az | (1:98\ | 0.22 |/0.97\| 0.27 | 0.41 | 0.20 | 
=) (on) da flvsas)] vsas | ative | ons | Sime@7+A tan 
2 | 0-19 | 0.1254 0.08 {0.08 [0.08 | emt bk 0 eee 
3a, v3 vs 3us | v3,u2 | w3,u3 sin (37+ A + 5x} 


YinGs/) Yar¥e | vide | b3,da | v2,u2 | va,y2 | Sim(T+2A +403 


Ys | YS | ves | ones | vege | vag | SiMOT+2A +508 
Pi] P? } pe fps P3 | Ps 
az | 0:32") 0.82 | 0.18") 0.18 | ==} 


Marte | Warhe | WalQe tebe ie Nee) cee sin(T+3A + a8) 


+ | MEY | obas | oder |] — | — | sin@T+3A +933) 


1Ay"P,° COs r+ 1B,"P,° sin v the 
2Ay'P;'cos (27+ A)+2B,'Py'sin(2T+ A), 
1Ay’P,* cos(T'-+2.A)+ 1B? Ps? sin( T+ 2A), 


1Ay Py’ cos T+ 1ByYP,°sin T, 
2A;'P3!cos (2 T+. A) +.B;'!P3'sin (2 r+ A), ( 7 ) 
2A;'P,%cos(2T+34 )+2B;°P,3sin (27434), 


are so conspicuous that we cannot neglect them. Figs. 2A and 2B show the variations 
of the pure universal time field. Figs. 2C to 2J illustrate the equipotentials of each 
i t : i iy vd Lt et 
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residual field at T=0, in which the component field }}{2.A,,'P,,!cos(2T+ A)+2B,'Pq' 
sin(2T+ A)} is most conspicuous. This is also expressed. as follows: 


dH2A,'cos(2T+24—4)+2B,'sin(2T+24—A)}P,! 
= 3124, cos (2#—A)+2B,,' sin (2#—A)}P,' 
=dha,!P,,'sin{2+20,(A)}, 
where : 


" F, . 2A, cos A—2B,,'sin A 
- 1\2 — = 1\2 rs 1\3 % = 22 n 
(s@, ) (A, ) +(.B, ») > tan; 2h, W, A= gf1. 1 sin A + Bs cos A > 


which may be regarded as a local time fiield of a half day varying with longitude. 
Fig. 3A shows the equipotentials of the total residual field at every successive two 
hours as compared with those of the local time field (Fig. 3B). 

Table 3 gives the 
corresponding values 


calculated by T. Naga- 


ta and M. Sugiura [7] B a Ve Nui ai ° a 
a ES V iaeN \\Y 
(and transformed to the {CANINES errt 
geomagnetic coordinate e V=Z Yea 
L YY SSZZzZ 


by the present author) 


on the basis of the 
dynamo-theory using 


SS tit 
peer, 
the result of Chapman’s THe ANS 


_ analysis of S;. Compar- 130° 


aye “ee 


ing Table 3 with Table Fig. 3B—Equipotentials of the local time field. 


2A, the main residual 

terms show a fairly good correspondence, so that the most conceivable cause of the 
residual terms seems to be an effect of the discrepancy between the geomagnetic and 
geographic axes upon the dynamo-action in the S,-layer, as shown by T. Nagata and 
M. Sugiura. In detail, however, this correspondence is not alway quantitatively 
sufficient, and it is not clear as to whether this variance is due to the assumption of 
a constant conductivity and of a velocity potential (using ¢,', ¢.? and ¢;°) in 
the calculation or to some other cause. Table 4 gives the contribution of the velocity 
potential to the magnetic one in the calculation based on the dynamo-theory, under 
the assumption of a constant conductivity. From this Table it seems likely that not 
only the velocity potentials ¢,', ¢,” and ¢;* but also ¢2!, 22, $43, %3' and $4? play a 
fair role in the S,-Field. Further detailed studies on the origin of the residual field 
will be reported in a subsequent paper, and we try, here, by calculation to obtain a 
coordinate system which is in good agreement with the distribution of the actual S,- 


variation. 
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3. A coordinate system for the S,-Field 

| By transforming the geomagnetic 
coordinates (9, A) to other suitable coordi- 
nates (6’, A’) having a new pole (9), Ao) 
(Fig. 4), the spherical surface harmonics are 


transformed as follows [9]: 
P (cos 9) cosmA = > —1)"[Aprt,»P%(cos 8’) 
w=0 
cos A! + Bun SnP% (cos 9’)sin 2A’), 


P (cos 9) sinmA=>\(—1)"[AypS,,P! (cos 6’) 
y=0 


Cos A! —Bym€mP% (cos @’)sin 2A’), ° 


Fig. 4 
where oD 
Ayn = Anp = a Gath (A +B), 
Bum = Bra = SV Oajden (A =}. 
d” EP Ac) 
_ —p\hem—-h* 1 “ 
A=(1—c)*s ar « +e) ee 
mp a” d*P,(c) ’ 
B=(1-+c)"s ae | al ——— ; 
P,,(c)=Legendre functions, 
Cmn=COSMAo, S,=SiNMA, 
C=cos@), s=sinQ, 
—m)! 
Gam = Ow eT e&o=1, p= Eg=*s+- Se, 2, 


By, this formula the general term in (1) is expressed as follows by referring to the 
new coordinate: 


tAnPn (cos @) cos (17+ mA) +,B™ n(cos@)sin(/T+mA) 

= 5 —1)"P" (cos O'(Aun—Bum) (Ac, +,Brs,,)cos((T+ nA’) 
+ (Aum—Bum (Bren —.Ams,,) sin (T+ pA’) 
+ (Aum t+Bum\cAney, +,B%s,,) cos (IT—A') 
+ (Aum +Bum\Bren—,As,,) sin IT — pA] 

=ELAwPs (cos 6’) cos (/T+-.4’)-+-,B™ Px (cos 6’) sin UT +- 2A’) 

+ CR’ P* (cos 0')cos (IT — pA’) + Dn Px (cos 0’) sin(IT—pA')] (4) 
where the terms ,C” and ~Dn in (1) were neglected. 

Since the new coefficients An’; Br’, .C™ and ,D™ in (4) are the functions of 
the new pole (05, Ay) respectively, we can illustrate the distribution of the sum of 
the squares of the 12-coefficients in the main Six residual terms on the 0) Ao-plane 
(Fig. 5). From Fig. 5 we find a minimum point (O,=<47, Ao = —170.8°). Hence, 
representing the S,-field by referring to the new coordinate system, which has this 
point and its antipodal point as the north and south poles respectively, the residual 
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MINIMUM POINT : (@,=4.1, A,=-170.8°) 


Fig. 5—Distribution of the sum cf the squares of the 12-coefficients in the main six residual terms. 


field become least as a whole. This coordinate system, therefore, seems to be more 
suitable for the S,-field than the geomagnetic one. 

The new coordinates (#’, A’) and meridian angle (¥’) of each station list the last 
three colums in Table 1, and as an example using this coordinate, we show a latitu- 
dinal distribution of the annual mean values of the range of 4H in Fig. 6. 

The north and south poles of our coodinate system are (82.5°N, 74.0°W) and 
(82.5°S, 106.0E) respectively (Figs. 7A and 7B), and it would be a question whether 
there is any significance that the north pole coincides with the centre of the auroral 
zone estimated by E.H. Vestine [10]. 

4. Conclusion 

From an analysis of the date during the Second Polar Year, we have © Remnie 
- a detailed feature of the S,-field, and found that this field contains a fairly conspicuous 
residual field in addition to the local time field. This residual field results in the 
variation of S, with longitude, and it seems that this field arise mainly from an effect 
’ of the obliquity of the earth’s magnetic axis with respect to its rotational axis upon 
the dynamo-action in the S,-layer. 

We have tried to derive a more suitable coordinate for the S,-field, by using a 
- coordinate transformation method to make the residuals least as a whole. Its north 
and south poles are (82.5°N, 74.0°W) and (82.5°S, 106.0°E) respectively. 
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An Average Equator for the Geomagnetie S,-Field* 


By Hireshi MAEDA 
(Geophysical Institute, Kyoty University) 


Abstract 


In the preceding paper [1] we obtained a coordinate system for 
the geomagnetic S.-field. The equator of this spherical coordinate 
may be regarded as an average equator for the S,-tield. 

This equator is in good agreement with a line of eastward 
maximum current density calculated by the dynamo-theory taking 
account of the discrepancy between the geomagnetic and geographic 
axes and of the high electric conductivity in the equatorial zone. 

These results seem to suggest the fairly deep-rooted effect of 
the dynamo-action due to the discrepancy of above two axes upon 
the S.fielkd on one hand, and the high electric conductivity in the 
equatorial zone on the other. 

1. Intreductien 

Since an anomalous daily variation of horizontal intensity at Huancayo had 
been pointed out by HE. Johnston and A.G. McNish (2] this interesting problem of 
diurnal variation of geomagnetism in equatorial regions has been discussed by many 
workers, M. Hasegawa has for a long time entertained the idea that these anomalous 
regions encircle the Earth, and he, with M. Ota, estimated such a line (say MS- 
egeefor) from an analysis of the data for the Polar Year, 1932-33 [3], Recently A.T. 
Price and GA. Wilkins [4) have also estimated a line of maximum 4H. These two 
lines are in good agreement within the limits of error (Fig. 1. Asa result of tem- 
perary observations in Africa, India and America, J. Egedal [5] has pointed out that 
these anomalous regions are situated along the magnetic equator. But from a result 
of temporary observations in south India and Ceylon by B.L. Gulatee (6], and by SK. 
Pramanik and SY. Narayanan {7} it should seem that the line of maximum 4H does 

If S, is due to an electric current flowing in the so-called Se-layer, there must 
be expected to be an anomalous concentration of electric current ( 
Chapman [8] in these regions. On the origin of an equatorial electrojet, A.G. McNish 
{9} attributed it to an effect of the discrepancy between the magnetic and geographic 
quators upon the dynamo-action in the S.-layer, and D.F. Martyn [10] to an increase 
of the electric conductivity near the magnetic equator. Recently M. Hirono [11], K. 
Maeda [12} T.G. Cowling [13], DF. Martyn and W.G. Baker {14} have discussed, 
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Say electrojet by S. — 


An Average Equator for the Geomagnetic S,-Field td 
independently, the conductivity of the ionosphere from a same point of view (Hall 
polarization mechanism), and given an explanation for the enhanced conductivity near 
the magnetic equator. 

We find here a line of eastward maximum current density 2t 11h 16m loczl 
time, at which H is nearly maximum, on the basis of the dymamo-theory, tzking 
account of the discrepancy between the geomagnetic and geographic axes and of high 
conductivity in the equatorial zone, and compared with the equator of our coordinate 
system for the S,-field. 

2. An Average equator for the S.-field 
The equator of our coordinate system obtaimed in the preceding paper [I] may be 
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regarded as an average equator for S,, 
3 Ip 30 190 and is illustrated in Fig. 1 as compared 


| 
: PEE eee bere with other several equators. This aver- 
Pa | a ea ae age equator for S, is in good agreement 
ff pe with that derived from the dynamo- 
Et Baal theory as below. 

$ The notation used is as follows: 
0, }=geographic colatitude and east 

longitude, 
r= geornasnetic local time, 
@, A=geomagnetic colatitude and east 


longitude, 


‘DEPARTURE PROM GEOGRAPHIC EQUATOR 


T=geomagnetic universal time (i.e. 


—69° meridian time), 
U, V=southward and eastward com- 
ponents of air velocity, 
¢=velocity potential of air, 
A=electric current function, 


“15° 
Fig. 5—Departure of the line of maximum cur- 
rent density @, of the equator of our 


coordinate system @ and of the geo- K=electric conductivity, 
magnetic equator (3) from the geographic H,=vertical component of the earth’s 
equator on the meridian of A =0. magnetic field 

> 


a=radius of the S,-layer. 


It is well known that the current function R for the stationary state must satisfy the 


following equation in the case of k=const.: 


OR 0 (6 9R\__,-{ (VH,) , 8(UH,sin@) 
“sindoae + a6 (sin 056) ec Las > a + ger s+ 


where H,=Ccos@ and C is a constant. If U and V can be derived from the velocity 
potential ¢, and ¢% is assumed to be expressed by surface harmonics. 


p= LW 


= 213) Rn PX (cos 0) sin (mr+ OS Sate in geographic coordinate, 
= pap Rnd) ch Pr (cos 9) sin (44+ Bx) eee in geomagnetic coordinate, 


(see our preceding paper [1], p. 48) 
where 


Cm= Ao, Sin(mT+-a™), 


. o 1 
Cn={ Alum Sin*(mT + an)-+ Bam Cos(mT +. as)" 
7 


2 , 
tan Bk = — fim tan(mT+-a™), 
Lum 
A,,, and B,,, =notations used in our preceding paper [1]. relation ¢). 


then the solution of eq. (1) is 
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Rs KCL dnd LemSn + Cir 1b + Crrq" sind + Bm), (2) 
where 
5 
St = n(n+-1y epee al ph (20-+-1)[2” + J 
(w+1) (2n+1) Phat n(2n+1) my 
sin*@ sin’@ 


p* = +cos 6)” q* = (cos 6)” 

Ct, C?;=arbitrary constants. 

Assuming that the electric conductivity 

of the layer concerned is fKy for region II 

and K, for regions I and III (Fig. 4), the 
solution for these three regions is 


Rr = —KoCLUD nen Sh +Ctp*)sin(uA @, 
4 +p), ---GEOMAG, EQUATOR-- --J[--- ~~~ 
Ru =— FRoCL Ln LenS +Cip*+ Crrrq*) gad 
sin (44+ 8m), 
Rin= —KoCU Dhan demSin +Cfyq"]sin(uA 
+ Br), 
(3) 


and using the condition that the normal Fig. 4 


component of the current and the tangential component of the electric field are con- — 
tinuous at the two boundaries of the three regions, the arbitrary constants are 


determined as follows: 


a Mu Au A. 
L aay x Cu= vi ’ Cin an 4a ’ Ciy = 4 ’ 
where 
A*=(f—-1" —M**(f4+ 1, 
=(f—1)c% We{(—1) #42 {M4 —(f-1)—-M*(F+1)}, 
ae, =(f-Dee Wey MF +1) —(F-D)1 
A =(f—Dek We LM (f+ D—(-D (Ff -D) 
At, =(f—1)c W% 2 FM —(— 1pm (fF — 1)+M**(f+1)}1 
~~ 1+ cos 9, x3 = S#(G1) 
Me-(itser) Be: 
I 
sew From the assumption of constant conductivity of the 
Ve=1on zen lan ok adaliig®y layer, the values of K,Ck can be estimated from 
. - | ma the observed magnetic potential [15]. Since the 
n nn n 
velocity potential ¢% contributes to the ‘current 
I 2 7.227 101° function -R%,;, a solution in region IJ, using the 
2 3 3.707 105° result obtained by M. Hasegawa and M. Ota [3] 
Zire 4 1.107 125° cabs 


(Table 1), becomes 
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(After M, Hasegawa and M, Ota) 
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x | 3p An Ban) sin (77+ “A+ Cn+t) 


= 5 (Anu+ Bun)sin(n Tp ter) | ? (4 ) 
where 

“ 1 927-E3 yf ae * 

st An n+2 () MoE n+ . 


Fig. 2 shows the lines of eastward maximum current density an 11h 16m local 
time in the cases of f=1, f=5, and f=100, and Fig. 3 the variation of current density 
along its line in the case of f=5. Fig. 5 shows a departure of the line of maximum 
current density from the geographic equator on the meridian of A4=0. From these 
Figures we have found that; 

(i) a line of the maximum current density for 0,=80° (as an example) and f=5 %s 
in good agreement with the above-mentioned average equator for S., 

(ii) the value for current density along the line of its maximum is large at the 
places where the geomagnetic equator is situated far from the geographic equa- 
tor, but the range of variation with longitude is merely 1.4 per cent, 

_ (iii) unless the highly conductive zone II is not very narrow, the position of maximum 
density may be shifted fairly well from the middle of the highly conductive zone. 

3. Conclusion 

Above considerations seem to suggest the fairly deep-rooted effect of the 
dynamo-action due to the discrepancy between the above two axes upon the S,-field 
on one hand, and the high conductivity in the equatorial zone (several times higher 
than in the middle latitude) on the other. 
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